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Anomalous spectral shift of transmission peaks is observed in a Fabry-Perot cavity filled with 
a chiral anisotropic medium. The effective refractive index value resides out of the interval 
between the ordinary and the extraordinary refractive indices. The spectral shift is explained 
by contribution of a geometric phase. The problem is solved analytically using the approximate 
Jones matrix method, numerically using the accurate Berreman method and geometrically using 
the generalized Mauguin-Poincare rolling cone method. The o-mode blue shift is measured for 
a 4-methoxybenzylidene-4’-n-butylaniline twisted-nematic layer inside the Fabry-Perot cavity. 

The twist is electrically induced due to the homeoplanar-twisted configuration transition in an 
ionic-surfactant-doped liquid crystal layer. Experimental evidence conhrms the validity of the 
theoretical model. 
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I. INTRODUCTION 

Optics of liquid crystals (LCs) is well known for fruit¬ 
fulness in applications and a remarkable variety of con¬ 
nections between observable physical phenomena [l[. A 
fascinating connection can be traced between the concept 
of geometric phase (GP) Q, also known as a topological 
phase, and a number of phenomena in quantum, relativis¬ 
tic, classical physics Q, and, in particular, in polarization 
optics p|-[^ Today photonics is at the apogee of topolog¬ 
ical ideas On the one hand, it originates from the 

condensed-matter graphene idea and the concept of topo¬ 
logical insulators [9j. On the other hand, it arises from 
the optical technology advance. For instance, recently 
in polarization optics the three-dimensional structure of 
the field of light polarizations with nontrivial topology 
has been directly measured 0- _ 

GP in twisted-nematic polarization phenomena m 
has application for the design of wave fronts using 
Pancharatnam-Berry phase optical elements [l^, [13 ■ 
Remarkably, polarizational GP is independent of the to¬ 
tal phase, so it is used to offset the frequency of a laser 
beam by GP modulators, adding the mechanical rota¬ 
tion frequency of a quarter-half-quarter-waveplate Pan- 
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charatnam device [14l - |l6l| . It permits switching by ferro¬ 
electric LG 0. 

Considerable attention is attracted to the research of 
LG placed inside a Fabry-Perot cavity (LC-FPC), com¬ 
bining small-voltage LG manipulation and high spectral 
resolution of the Fabry-Perot interferometer. Fundamen¬ 
tal photonic degrees of freedom—in transmittance [H- 
[2l|, phase [13, and polarization 231—can be efficiently 
controlled. Polarization control usually uses the Mau- 
guin adiabatic waveguide mode (23 | in a chiral anisotropic 
medium, particularly in a twisted-nematic (TN) LG layer 
inside the Fabry-P&ot cavity (TN-FPC) [l^ . 

The basic method for TN-FPC calculation is the Jones 
formalism of complex vectors and matrices of dimen¬ 
sion 2 [ 13 . Abele [13, HI] had introduced Chebyshev 
identity for the matrix power which was successfully em¬ 
ployed both in layered and in anisotropic media [13- 
[3l|. To solve the problem one finds an eigenwave (op¬ 
tical mode) which conserves its shape while propagating 
through the medium. A set of eigenwaves is described by 
the eigenvectors and the eigenvalues of the Jones matrix. 

Another approach to find the eigenwaves is to solve 
the Riccati-like ordinary differential equation system [sH - 
113 • Within the framework of coupled-mode theory and 
modal analysis, this approach is equivalent to the matrix 
one (see Appendix C in [13 )■ 

The account of anisotropic reflections between the sub¬ 
layers of bulk LG led to a generalized Berreman formal- 
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ism for matrices of dimension 4 [Sdl - fs^ . This generaliza¬ 
tion is necessary in media such as a cholesteric LC [39l - l4^ 
and a TN cell of small thickness, and other media with 
a sharp spatial variation of dielectric characteristics at 
the wavelength scale However, Jones formalism 

gives a good approximation when the thickness of the TN 
cell is several times larger than the wavelength and the 
dielectric constant varies smoothly. 

Assuming no reflection between sublayers of the LC 
bulk, the TN-FPC behavior was described at high volt¬ 
ages [ 2 ^ as well as at low voltages (S^, ■ The connec¬ 

tion between these two extreme cases is described in 
and generalized in [d^. Another approach is the substi¬ 
tution of the multilayer medium by an effective homo¬ 
geneous anisotropic pla te [5 ll| . The independent method 
is incorporated in [24l . l52l| using distinct mathematical 
tools: group theory and phase space. Also TN-FPC can 
be considered as a one-dimensional photonic crystal [H- 
[ 5 ^ . The photonic crystal itself can be formed by a LC 
material 

The orientation model of the LC layer is required to 
determine its optical response. To the best of our knowl¬ 
edge, in a TN under electrical voltage the general orien¬ 
tation model cannot be derived analytically and is simu¬ 
lated numerically, even for the one-dimensional case . 
In contrast to in-plane-oriented nematic [^ . TN pro¬ 
duces optical mode coupling manifested as avoided cross¬ 
ing of spectral transmission peaks . A certain way to 
eliminate the mode coupling using anisotropic mirrors is 
suggested in M- The original theoretical study of the 
apparent paradox of the mode number j ump for mode 
coupling inside a TN-FPC is proposed in (4^ . 

This article examines the TN-FPC sample with no 
electric voltage deformation that allows analytical de¬ 
scription. The direction of the spectral shift of the trans¬ 
mission peaks is far from obvious. That’s why a visual 
connection is presented for the proposed spectral shift 
and the geometric phase shift. This shift observation is 
hindered by interplay of four optical waves of opposite 
directions and orthogonal polarizations. Positive feed¬ 
back condition describes the total spectral shift (SS), 
assuming three types of wave coupling. The first type 
of the coupling originates from LC twist which induces 
a twist spectral shift (TSS). The second coupling type 
is produced by cavity mirror reflection which induces a 
reflection spectral shift (RSS). The third coupling type 
is produced by reflection between sublayers of LC bulk; 
in this research it is assumed insignificant. Previously 
developed theory [s^, is generalized for the account 
of anisotropy at mirror reflection, meaning the distinct 
reflection phases of e- and o-wave components. The re¬ 
flection is anisotropic even when the mirror is made of 
an isotropic material while the cavity medium itself is 
anisotropic. 

The experimental scheme excludes significant spectral 
shift impact from parasitic factors other than TSS and 
RSS. The experiment confirms the theory qualitatively 
and quantitatively. 


II. ANALYTIC MODEL 


twisted 

nematic ^^"'•'0'' 


FIG. 1. Sketch of TN-FPC, a cavity with a chiral anisotropic 
medium. 

A FPC consists of two plane mirrors (Fig. [T]). The 
reflecting surfaces are oriented in xjz-plane. Nematic is 
placed between the cavity mirrors. The LC director is 
a unit vector of predominant direction of LC molecules. 
Twist is the state when nematic layer is divided into thin 
lamellar sublayers with the nematic director being con¬ 
stant inside every sublayer and rotating from sublayer to 
sublayer. Uniform twist with no orienting external fields 
is implied when the LC director rotates uniformly in the 
plane of sublayers along the right screw. In Fig. [T] the 
twist angle is 80 degrees, hence the analysis is valid for 
an arbitrary angle. The LC director field determines the 
local dielectric tensor all over the medium. The extraor¬ 
dinary dielectric permittivity axis is collinear to the LC 
director. Consider the nematic with a positive uniaxial 
anisotropy. The extraordinary and ordinary refractive in¬ 
dices (RI) correspond to waves with slow and fast phase 
velocities and equal to ne,o = n± 6n. 

Let the average phase cr, anisotropy phase (angle) d 
and twist angle tp be linear functions of the coordination 
z along the layer normal direction: 

a (z) = nko z, S (z) = 6k z, ip (z) = k^ z, 

where fco = w/c is the angular wavenumber, 5k = 6nko, 
kip = ip (L) /L, L is the nematic layer thickness or dis¬ 
tance between mirrors. At z = L let function values 
a {z = L), 5 [z = L), ip{z = L) be written simply as cr, 
6, ip without the function argument. Let the light im¬ 
pinge to the TN-FPC strictly in z-direction. Electric field 
strength is described by the pair of x- and ^-projections: 

Ex (z) exp [i {ujt — a (z))] -I- c.c., 

Ey (z) exp [i {(jjt — a (z))] -I- c.c., (I) 

where c.c. stands for the complex conjugate component. 
The pair of complex amplitudes E^^y are convenient to 
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be written as the Jones vector 


&xy (z) 


Ex {z) 

Ey (z) 


The Jones matrix for untwisted nematic at </? = 0 is diag¬ 
onal when the nematic director is collinear to the a;-axis: 


A(J) 


0 

0 


( 2 ) 


Note that the extraordinary field component along the x- 
axis has negative phase shift according to standard notion 

O: 


permittivity main axes of neighboring sublayers. Itali¬ 
cized form “d” is used to distinguish it from the particu¬ 
lar case of infinitesimal differential operator. 

The rotation matrix is written as 

cos(p sinyi 

— sin (f cos (f ' ' ' 

It rotates the reference frame about the z-axis by the 
(f angle. The Jones matrix for the waveplate (retarder) 
situated at (p angle is written as: 

Ay, = 

In rotating frame the polarization ellipse appears to be 
rotated by negative angle —(p. The total Jones matrix is 
written as the matrix product: 




^xy (E) — ACxy ( 0 ) 


Ex ( 0 ) 

Ey i0)e+^^ 


The Jones matrix A is the transfer matrix (or propaga¬ 
tion matrix) of the layer. It transfers a polarization state 
from one layer boundary to another by multiplying the 
corresponding Jones vector. 


A. Traveling eigenwaves for a uniformly twisted 
nematic 


Jo = $ 

X $ 



A{d6)... 


(^)aw 4 



X 4 A (dS) 4 = 4 (-f,) J. 

The product is supposed to be read from right to left 
with increase of z, because the Jones column vector is 
substituted on the right side of the matrix. In rotating 
frame of the matrix J the LC director is always collinear 
to the primary axis. It is the e — o-frame which is often 
used in description of chiral media [s^, . 

z-axis rotations are additive, 


For convenience of subsequent interpretation let us 
present a direct trigonometric derivation of some general 
expressions for a traveling wave in a chiral anisotropic 
medium. These classical results can be verified in [3 If . 
Let us divide a TN layer into a series of equal sublay¬ 
ers each with the thickness dz and the anisotropy an¬ 
gle dS = 6 (dz) = 6k dz. Generally, the following is 
valid for finite not twisted sublayers. The twist angle 
dp = ip {dz) = kip dz is the angle between the dielectric 

I 


l> {p2) ^ {pi) = ^{p2 + Pi), 


so the Jones matrix is naturally decomposed into the 
product of certain sublayer matrices: 


J = d = 


dp 


$(^)A(dJ)cI>(^l 


dp\ 


Ns 


( 4 ) 


where Ns = p/dp is the number of sublayers. Substitu¬ 
tion of Eqs. @ and m gives: 


cos{dp/2) sin(d(^/2) 
— s\n{dp/2) cos {dp/2) 


^-idS Q 

0 e- 


cos {dp/2) sin {dp/2) 

— sin {dp/2) cos {dp/2) ’ 


dJ = 


cos {dp) cos {dS) — i sin {d6) sin {dp) cos {d6) 
— sin {dp) cos {d6) cos {dp) cos {d6) 


i sin {dS) 


Eigenvectors ej of the matrix dJ describe polarization 
conserved in the rotating basis, and eigenvalues gj of the 
matrix dJ are phase factors of propagation through the 
medium layer. The eigenvalue condition is: 


det 


= 0 , 


(dJ - gji^ 
det {dJ^ — tr {dJ^ gj + gj = 0, 


( 5 ) 


( 6 ) 


dJ ej = gjej 


where / is the identity matrix. According to Eq. ((S]) the 
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determinant of the matrix equals 
det (dJ^ = 1- 

Transfer matrix is unimodular when the transferred en¬ 
ergy is conserved. The trace of the matrix is 

tr (dJ^ = 2 cos {dip) cos {d6). 

Let’s introduce a new twisted anisotropy angle defined as 
the following: 

cos {dv) = cos {dip) cos {dS). (7) 

The solution of Eq. ([6]) can then be written as 
= cos {dv) =F i sin {dv) = exp {^idv). 

The eigenvectors take the forms: 


angle or the twist angle. With zero twist dip = 0 and 

= 0, the eigenwave ej is simplified into extraordi¬ 
nary e-wave, and ej is simplified into ordinary o-wave. 
They are named as quasi-e- and quasi-o-wa.ves |34l |. or te- 
and to-waves (“twisted waves”) . The latter italicized 

form is appropriate here. However elliptically polarized 
to-wave may be confused with linearly polarized TE-wave 
(transverse electric mode), whose electric field is perpen¬ 
dicular to the reference plane or axis. Also this abbrevia¬ 
tion (extraordinary/ordinary) may be confused with the 
parity abbreviation (even/odd). The chirality of to-wave 
is opposite to LC director chirality. The chirality of to- 
wave is the same as LC director chirality. Consequently, 
the wave pair use to be termed as “opposite chirality” 
wave and “same chirality” wave [b^. In cholesteric LC 
to-wave demonstrates the Bragg reflection, while to-wave 
is non-diffractive. That is another way to distinguish the 
eigenwaves 1^. 

The matrix dJ is diagonalizable using the unitary ma¬ 
trix U to transform basic vectors te^o into ete,to- 


Cj — Cie — 

cos'd 
—i sini? 

; Cj — 6to — 

1 1 

O) g 
u 

1 _! 

, (8) 

1 

II 

re 

p 

II 

cos'd 
—i sini? 

—isini? 
cos'd 

where 





^ ^ L ^to J 

COSl? 

i sinr? 

zsini? 

cosi9 


= 0 / 2 , 

cos 0 = sin {d6) / sin {dv ), 

sin 0 = sin {dip) cos {d6) / sin {dv) 


(9) 


In the literature the pair of eire nwaves described by (|S]) 
has several names [^, 1^, . As usual the terminol¬ 

ogy is chosen depending on domination of the anizotropy 


where the symbol “f” indicates the Hermitian conjuga¬ 
tion. 


dJ = U- 


-idv 


0 

„-i-idv 


u. 


( 10 ) 


The diagonalization simplifies the Jones matrix (Eq. (|3])) 
exponentiation: 


J = 




Ns 


= <u 


y-1 


— idv 


^-\-idv 


Ns 

=U-^ 


— idv 


n Ns 


a+idv 


u = u 


-1 


a-^iv 


U, 


J = 


cos'd — fsin-d 

-fsin-d cos'd 



■ g-™ 

0 



0 

g+it- 



J = 


cos V — i sin v cos 0 sin 0 sin v 
— sin 0 sin v cos v i sin v cos 0 


■ (11) 


The uniform twist condition was used: 

Ns = ^ = - = — = -. 
dip dS dv dz 

The solution (HID is valid for sublayers of finite thick¬ 
ness and is the representation of Chebyshev identity [2^ . 
Equation 0 for the angle dv is the Pythagorean theo¬ 
rem for the spherical right triangle. Eor the smooth twist 
function ip {z) the sublayer thickness dz tends to vanish. 
The solution can be simplified using the Pythagorean the- 


cos d f sin 
i sin 'd cos 'd 


I- 

orem for the plane right triangle: 

dv^ = d5^ -\- dip^. 

It is easy to derive it by tailoring the cosines in Eq. 0. 
Multiplication by Ng produces 


2 r2 , 2 

V =0 ip ■ 

The total wave phase then is written as 


( 12 ) 


CT ± r; = cr ± sjdP- ^ ip^. (13) 

Let’s name it as the Mauguin formula. The effective RI 
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is found by dividing both sides by k^L: 

nte,to =n±\J S-n? + [ip/k^hf. (14) 

The eigenwave ellipticity parameter 0 from Eq. can 
be reduces as relation 

tan© = (f/S. (15) 

Physically this ellipticity is the smoothness of the twist 
angle growth in comparison with the anisotropy angle 
growth. It is the adiabatic parameter of Mauguin’s 
waveguide regime. 


Assume the half-turn rotation matrix to act on the axis 
perpendicular to LC director. 


R = 


-1 0 
0 1 


Obviously, the double reflection of produces the iden¬ 
tity matrix. 


C. Perfect cavity eigenwave 


B. Mirror reflection matrix 


The mirror reflection matrix has two multipliers: the 
phase marix Mq and the half-turn rotation matrix R. 
For a mirror made of metal with R1 rim the reflection 
originates from high R1 contrast with LC {rim — ri ^ 
Sn). In this case the phase matrix Mq is approximately 
isotropic. However for the dielectric multilayer Bragg 
mirror, the R1 contrast is not high and the phase matrix 
is far from isotropic. 


Mn = 


exp(-i/ie) 0 
0 exp {-ifJ.o) 


= exp (-*CT^) 


exp(-zd^) 0 
0 exp i+iSfj,) 


= exp {-iUfi) M. 


Various algorithms are used to find phases /ie,o for certain 
mirrors 


Assume the cavity is perfect (free of losses). The whole 
loop of the wave propagation through the cavity consists 
of a couple of passages and a couple of reflections. The 
corresponding matrix is the square of the half-loop ma¬ 
trix: 


L = exp {—i(7 — 

= exp {—2ia — 2icr^), (16) 


where H = RMJ. 

If the polarization matrix H is not an identity matrix, 
then its eigenvectors coincide with the eigenvectors of L. 
The eigenvalues of L are expressed through the eigenval¬ 
ues of H: 


9l = gn exp {-2ia - 2icr^). 


Now one can find the eigenvalues of H: 


H = RMJ = 


J 


■-1 o' 

exp {-iS^) 

0 

cos V — i sin v cos 0 sin 0 sin v 

0 1 

0 

exp (-1-1(5^) 

— sin 0 sin v cos v + i sin v cos 0 


H = 


— exp {—iSfj,) (cos t> — i sin u cos 0) — exp {—i6^) (sin 0 sinu) 

— exp {+iS^) (sin 0 sinu) exp {+iS^) (cos v + i sinu cos 0) 


det = —1, 

tr = 2i (cos u sin (5^ -|- sin u cos 0 cos d^) = 2i cospo- 
The eigenvalues are 

gH = i (cospo T * sinpo) = exp {i (7r/2 =F Po)) = ± exp {±i, 
where the resonator phase p = 7r/2 — po, 

sin p = cos Po = cos v sin -|- sin v cos 0 cos (17) 


The solution corresponds to a couple of modes resonating 
inside the cavity (standing waves) and let’s denote them 
as re- and ro- waves. 

In the exotic case where 5^ = 7r/2 Eq. (flTl) is simplified 
to Po = ±u. In [6lj the anisotropic mirrors are suggested 
with phases pe = tt, po = 0; consequently, <Jfi = 5fi = 
7r/2. 

gL = exp{-2{a -b tr^ ± p)) 

= exp{-2{a =F Po) + ’’"(I ± 1)) 

= exp{—2{a ± v)). 

The reflection coupling is removed and RSS is eliminated. 
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In the isotropic reflection case of (5^ = 0 Eq. (El) 
reduces to 

sinp = sinucos©. (18) 

It is equivalent to expressions presented in [13, . 


D. Positive feedback condition 

The position of transmission peaks is determined by 
eigenfrequencies of perfect cavity waves (modes). They 
satisfy the positive feedback condition for the total 
phase shift to be a multiple of 27r. 

—2(7 — 2(7^ 2p = —2ttN, 

= a + - ttN, 

where N is the cavity mode number. Using Eq. ED: 

=F sin (t 7 + /i — ttN) = cos v sin + sin v cos 0 cos 

(19) 

Here the minus sign corresponds to te-wave, whereas the 
plus sign corresponds to to-wave. It is possible to solve 
this trigonometric equation graphically. 


E. Dispersion curves and TSS 


F. Dispersion curves and RSS 

Eigure[2](b) illustrates dispersion curves for re- and ro- 
waves under simplification by Eq. m- 

=F arcsin(sinucos 0 ) = (7 — ttA^. ( 21 ) 

The left hand side of Eq. (I^Tl) at cr^ = 0 was shown ear¬ 
lier in [131 as a resonance diagram. The right hand site 
of Eq. (1^ produces constant-slope lines. The reduced 
zone with the period n corresponds to the half-loop. The 
positive feedback condition m is fulfilled at every inter¬ 
section of the magenta curve with the green or blue one 
for re- and ro-wave frequencies corresponding to spectral 
transmission peaks. Dashed dispersion curves correspond 
to te- and to-waves =F arcsin (sinu cos (0 = 0 )) = ^v. 
These lines describe the zeroth-order approximation, no 
reflection coupling. Points T and R indicate the RSS of 
mode 3 from to- to ro-wave. 

Dispersion curves for te-, to-, re- and ro-waves meet at 
points Gi_ 2 ... where the Gooch-Terry minimum condition 
is valid: 

sin (v) = 0 . ( 22 ) 

The condition was obtained for minimal TN cell trans¬ 
mittance [33| . In Gooch-Terry minima both resonator 
phase p and twisted anisotropy phase v are multiples of 
TT. The Jones matrix J (Eq. (fTTl) ') is degenerated into a 
unit matrix: 


Eigure[2Ka) shows dispersion curves of te- and to-waves 
for if = 7 r /2 and ip = Q. Scales are dimensionless for 
both axes. At the ordinate the o-mode number No = 
{a — S) 121^ = 2L/Xo is proportional to the frequency of 
the light field. At the abscissa the phase shift is the wave 
number multiplied by the cavity length. 

The branch of to-wave does not show the splitting by 
the cholesteric stopband ( 0 , at p. 354), for tp = 7 r /2 
at t7 = u, in the point B of Fig. [D(a). The splitting is 
omitted in Eq. ED- 

Dashed curves for the untwisted structure correspond 
to o- and e-waves: v {p = Q) = 6. Points O and T indi¬ 
cate the TSS of mode 3 from o- to to-wave. Calculated 
curves for te- and to-waves lie outside of the intervals of 
o- and e-wave phases. It illustrates the fact that effective 
RI lies outside of the interval determined by the ordinary 
and extraordinary RI. 

Let’s use Eq. ED to determine the twist phase shift: 


j = ±i 

with ±1 eigenvalues and arbitrary eigenvectors. Let the 
condition 

sin^ V = 1 (23) 

be the Gooch-Terry maximum condition. Note that this 
simple condition describes local maxima of the TN cell 
transmittance only approximately. In fact, the Gooch- 
Terry transmittance is a sine function of phase, and max¬ 
ima of this function are slightly non-equidistant. Figure 
HDb) shows that in Gooch-Terry maxima the phase p is 
maximally distant from the phase ±v. This difference is 
given by substituting the condition (I23p into Eq. m 

V = arcsin(l) = 7r/2, 
p = arcsin (cos 0) = 7r/2 — 0, 

V — p = Q. 


V — S = \/(52 _|_ (^2 _ 
assuming <C J gives 



The maximum difference is described by the adiabatic 
parameter 0. Near the Gooch-Terry maximum the dif¬ 
ference decreases according to the law: 


min {tt — p — v,v — p) 


Vl J- 0^ tan^ V — 1 

|tanu| 


(24) 


V — S = 6 


( 20 ) 


Near the Gooch-Terry maximal point G 3 / 2 , the four 
dispersion curves form a typical pattern called avoided 
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FIG. 2. The dispersion curves. Abscissa axis shows the phase shift proportional to the wave vector. Ordinate axis shows the 
number of modes, proportional to the frequency of the wave. Blue color shows o-, to-, ro-waves; green color shows e-, to-, 
re-waves. The calculation parameters are ip = 7 r/ 2 , 5n/n = 0.3. (a) TSS calculated with Eq. (I13II . Points O and T indicate 
the frequencies of mode 3 for o- and to-waves, respectively. Point B corresponds to the zeroth mode. The splitting at the 
point B is not shown, (b) RSS calculated with Eq. (I21II . Points T and R show the frequencies of mode 3 for to- and ro-waves, 
respectively. Gi and G 2 are Gooch-Terry minima (I22II and G 3/2 is the Gooch-Terry maximum (I23II . 


crossing. In oscillation theory re- and ro-waves corre¬ 
spond to normal frequencies, while te- and to-waves cor¬ 
respond to partial frequencies [d^. In quantum mechan¬ 
ics re- and ro-waves correspond to adiabatic states while 
te- and to-waves correspond to diabatic states [t^ • Both 
terminologies are used in optics. 

The alternation of Gooch-Terry minima and maxima 
in Fig.[2Jb) can be interpreted as an alternation of cross¬ 
ings and avoided crossings of transmission peaks [13, ■ 

Conditions (1221) and (1221) for the untwisted medium cor¬ 
respond to the respective conditions of sin^ <5 = 0 and 
sin^ d = 1. The Gooch-Terry minimum matches the 
twisted analogue of a waveplate. This waveplate retardes 
te-wave compared to to-wave by an integer number Ns of 
wavelengths with phase shift 2t:Ns. Gooch-Terry max¬ 
ima correspond to phase shifts of {2Ns -\-1 ) tt. 


G. Spectral shifts 


And substituting p from Eq. m gives: 

AA = — --- (arcsinlsinucos©) — d). (25) 

2'KnoI^ 

Using approximations (|221) and dini) one obtains: 


AA = AAtss + AArss 


A^ 

2'KnoL 



Vl -1- 02 tan^ V — 1 


tanu 


(26) 


Far from the Gooch-Terry maximum condition, the sec¬ 
ond summand AArss can be neglected: 


AAtss = — 


AV' 


A 

2TrnoL \ 2S 
A3 


2Tmo2L 2'K6nL/\ 2ndn 




Consider the total spectral shift AA of the TN com¬ 
pared with an untwisted counterpart. Without loss of 
generality we consider the o-wave. For the e-wave one 
can derive symmetric formulae with the opposite sigh: 
AAe = —AAo 

Dimensionless relative shift of the vacuum wavelength 
A is: 

AA Ak p — 5 
A ko — Sk a — S 

The average phase a can be excluded by use of the fol¬ 
lowing relation: 

A _ A2 
a — S 2TTnoL 


For the o-wave TSS is wavelength-negative AAtss < 
0. And peaks are shifted in shortwave region; i.e., 
blue-shifted. However near the Gooch-Terry maximum 
tan (v) —>■ 00 so that 

^ ip^ \/l -I- 02 tan^ V — 1 pP' ^ (p^ — 2ip 

2S Itanoj 2S ~ 2S 

For example, at (^ = 7r/2, 

TT — 4 

AX^p-2= -< 0. 

2 

Thus at p < 2 (in radians) the RSS component may 
locally reverse the total shift direction from blue to red. 























(a) 


(b) 



FIG. 3. Representative polarizations at medium boundaries z = 0 (a) and z = L (b). The twist angle ip is between o- and 
^/-directions. The o-defiection angle ^ is between o- and ro-directions. The e-deflection angle between e- and re-directions has 
the same magnitude as ^ at both boundaries. 


TABLE I. Three levels of anisotropy complexity. 


Anisotropic medium Homogeneous Twisted 

TN-FPC 

Eigenwave 

o,e to,te 

ro,re 

Eigen 

Linear Elliptic 

Linear 

polarization 

at boundaries 


III. INTERPRETATION 

Three levels of anisotropy complexity are presented in 
Table H] and illustrated in Fig. [3] Eigenwaves re- and ro- 
have linear polarizations at mirrors of the cavity [U13 ■ 
These polarizations are biased from the LC director and 
the direction orthogonal to the LC director by the deflec¬ 
tion angle With increasing frequency linear polariza¬ 
tions rotate continuously passing from one of the princi¬ 
pal axes to another. Therefore, we suggest to denote the 
eigenwaves as e2o- and o2e-waves (English words “two” 
and “to” are pronounced identically). At the Gooch- 
Terry maximum the polarization coincides with the bi¬ 
sectors ^ = ±45°. Therefore, the waves are called bisec¬ 
tor and orthogonal bisector [3 • It is convenient to treat 
e2o-wave as re-wave while it is close to e-wave and to 
change its name to ro-wave after it passes through the 
bisector and gets close to o-wave. And vice versa for o2e- 
wave. Dispersion curves in Eig. [2](b) show this renaming 
by changing color at the Gooch-Terry maximum. The 
curve G 1 RG 2 for o2e-wave is blue-colored in the lower 
part for ro-wave. And it is green-colored in the upper 
part for re-wave. 


A. Poincare sphere 

There is a variety of methods to make the results more 
comprehensive and visually attractive. They are the 
complex-number representation of polarization, Poincare 


sphere [3l|, [H, [3, high-order and hybrid-order 

Poincare sphere [73|, [3, admittance diagram method, 
Volpert-Smith chart and three-dimensional (3D) Smith 
chart [iilzl, and the rolling cone method [2J, |3 ■ The 
last method provides a mechanical visualization of Eqs. 
(HU and © with a sum of orthogonal angular velocities 
of a solid cone rolling on a plane. Por the sublayers of 
finite thickness the cone is replaced by the pyramid. The 
cone is rolling in a characteristic space of light polariza¬ 
tion ellipses. This space is usually represented as a sphere 
of unit radius called the Poincare sphere (PS). 

Eigure U shows a to-wave polarization trajectory and 
a ro-wave polarization trajectory (in blue and magenta, 
respectively). They smoothly evolve with the penetra¬ 
tion depth of the TN layer under the action of operator 
Jq. For the right-handed to- and ro-waves it is conve¬ 
nient to set the right-handed polarization in the upper 
hemisphere of PS as in Issj ]. and not in the lower 
one, as in [1, l23- The trajectory RR' is the spher¬ 

ical trochoid [3- It is associated with a trajectory of 
a point rigidly connected with a solid cone (Fig. [S](a)) 
rolling without slipping on a plane [3 - Stereographic 
projection of a similar trajectory is presented in Ref |33l |. 
p. 136, Fig. 2.24. 

For further consideration it is essential that the uni- 
modular Jones matrix has the geometrical sense of the 
PS rotation. Jones vectors correspond to PS points. Dis¬ 
tances between points are conserved at transformations 
on the spherical surface. The PS point is often regarded 
either as the normalized triplet of Stokes parameters or 
as the polarization ellipse traced by the terminal point of 
the field vector. In both cases PS is a two-dimensional 
manifold. But the Jones vector has three degrees of free¬ 
dom. Its third phase is the temporal phase. It progresses 
by 27r while the field vector passes the elliptic trajectory. 
The 3-phase polarization state is the one-to-one repre¬ 
sentation of the Jones vector. This polarization state is 
regarded either as the triplet of Euler angles or as the 
unit quaternion. The space of all polarization states is a 
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(a) 


(b) 




O, deg 


FIG. 4. Poincare sphere (a) and a portion of its cylindrical 
projection (b). Directions y, x — y, x, x + y correspond to 
angles <E> = 2(^=0, 90, 180 and 270 degrees, respectively. Tra¬ 
jectories under the action of operator Jo for ip = 80°: OO' 
- linear polarizations on the PS “equator” corresponding to 
o—wave, TT' - to-wave trajectory on the PS “parallel” with 
latitude 0, RR' - ro-wave trajectory, a spherical trochoid. 
The parameters correspond to the rightmost peak in the men¬ 
tioned below experimental spectrum of Fig. Eta), A = 579.1 
nm. 


3-sphere. The Hopf fibration projects it onto PS which is 
a 2-sphere The unit quaternion may be imagined as 
a “flag” consisting of two arrows. The first polarization 
arrow goes from the center of PS to its surface. The sec¬ 
ond arrow of temporal phase is connected to the terminal 
point of the first arrow and goes in a perpendicular direc¬ 
tion (similar to that shown in Fig. 2 of Ref [Z^). If the 
polarization arrow is rotated then the entire “flag” of po¬ 
larization state is rotated about the same axis. With the 
temporal phase increase the second arrow rotates around 
the first one. Remarkably, one period on the polariza¬ 
tion ellipse corresponds to two revolutions of the “flag.” 
Strictly speaking, Jones matrices form the special unitary 
group SU([2). This group is the universal covering of the 
rotation group SO ([3]). The covering is two-sheeted [sol - 

m. 


B. Geometric phase 

The parallel transport of a geometric object on a 
curved surface rotates the object about its own axis. A 
classic example is the Foucault pendulum with the rota¬ 
tion of its swing plane caused by the Earth’s daily ro¬ 
tation. Similarly, the parallel transport of the polariza¬ 
tion state on the PS curved surface leads to the phase 


shift called the geometric phase (GP). It is caused by 
global geometric characteristics, such as the curvature 
and the parallel transport trajectory, and independent 
of local characteristics, such as the speed of state move¬ 
ment along the trajectory. There is a geometric formula 
for closed trajectories. Applied to polarization optics it 
claims: GP /3 is equal to minus half the area ft encircled 
by the trajectory on PS: 

P = -nf2. (28) 

A rigorous proof is given in @ using Stokes’ theorem. 


C. Geometric calculation of the phase shift 
corresponding to TSS for traveling wave 


We apply the geometric formula to find the phase shift 
of the traveling wave at the trajectory TT' in Fig. 2] 
Elliptical wave expressed by Eq. m is a superposition 
of o— and e—waves with relevant RIs. The averaged RI 
should be chosen in the form that takes into account the 
analogue of the Aharonov-Anandan dynamic phase (see 
Ref 0, Eq. (5.19)): 


n = n —Sn cos Q, a = nkoL = a — 6 cosQ. (29) 

The minus sign corresponds to o-wave. This RI nor¬ 
malization reduces the PS rotation to parallel transport 
along the great-circle trajectory. 

Let the twist angle ip = n. Then <I> = 2(^ = 27r makes 
one turn on the PS “parallel” with latitude 0. The area 
between the “equator” and this “parallel” is equal to the 
side surface of a cylinder of unit radius and of height 
equal to sin© (see Ref [S^, Lambert cylindrical equal- 
area projection on p. 260). 

ft {ip = tt) = 2tt sin 0 . 

For arbitrary twist angle 

fl ((p) = 2(/3sin0. (30) 

The total phase shift 7 consists of the dynamic phase a 
and the geometric phase /3: 7 = a -I- /3. The geometric 
formula (E51) and equations (I29I30I) give: 


7 = cr — J cos Q — ip sin 0 . 

Eqs. (fT^ and (fT^ produce a transformation 5 = v cos 0, 
ip = t!sin0. This transformation gives 

'y = a — V (cos^ 0 + sin^ 0) = cr — \/ 5'^ + ip'^. 

Indeed, the result coincides with the Mauguin formula 

m- 
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D. Geometric calculation of the phase shift 
corresponding to RSS for cavity wave 

The Mauguin-Poincare rolling cone method is easily 
expanded to the Hopf bundle of PS. In fact, the gen¬ 
eralization is a consequence of the solid cone analogy. 
This generalization allows to find certain phase relations 
trigonometrically. The phase integration along the tra¬ 
jectory RR' corresponds to the matrix multiplication of 
Eq. (jT]) . Diagonalization of Eq. (flUl) is a simplifying al¬ 
gebraic transition to the rotating frame written as the 
product of three matrices. Geometrically, it corresponds 
to the transition from the trajectory RR' to the chain of 
three arcs: RT-TT'-T'R' (Fig. |Tl[b)). Wherein the first 
and the last arcs are to be chosen as geodesics. The great 
circle arc RT (see Fig. [IKa)) matches the parallel trans¬ 
port and the geometric phase corresponds to the Pan- 
charatnam phase 0 ] without any normalizing dynamic 
phase of Eq. (OH)) . 

During the passage of the distance L ahead through the 
cavity the traveling to-wave in rotating frame receives the 
phase shift v. Therefore, PS is rotated by the angle 2v. 
After the rotation the linear polarization R is returned 
back to the “equator” at the point R'. It corresponds to 
reflection symmetry of triangles OTR and O'T'R'. The 
mirror reflection corresponds to the triangle reflection 
symmetry with respect to the RO arc (improper rota¬ 
tion). The resultant geometric phase shift of the ro-wave 
compared with that of the to-wave accounts the triangle 
area four times per loop cycle 4D {OTR) = 2/3rss- 

The rolling cone method works even when the adia¬ 
batic condition is violated. In the non-uniform LC 
twist is examined, for example, under electric voltage. 
However, in the symmetric case for the cone opening 
angle 0 (z) = 0 (L — z) the eigenwave preserves linear 
polarization on the layer boundary. For practice the lin¬ 
ear polarization of the transmitted light is advan tag eous, 
because of the reliable blocking of transmission |^ . 

The account of mirror phase shift 6^ by anisotropic re¬ 
flection is equivalent to a waveplate action. On PS this 
leads to an additional rotation of the RO arc by the an¬ 
gle of Sfj, off the “equator.” Pythagorean theorem (fl^ 
for the right spherical triangle OTR transforms to the 
cosine theorem (El) for the angles of the spherical trian¬ 
gle. The point R leaves the PS “equator”. It means that 
the corresponding transmission peak becomes elliptically 
polarized. The superposition of opposing traveling waves 
on a perfect cavity mirror remains linearly polarized at 
yf 0 only when cr^ = ttN^/2 with integer N^. 


E. Connection between dispersion curve phases 
and PS angles 

With increasing frequency the Mauguin cone rolling 
angle 2v is increased uniformly. The to-wave corresponds 
to the Mauguin cone axis on PS. The point T is fixed, 
assuming constant adiabatic parameter: 0 (w) ~ const. 




FIG. 5. The rolling cone method on the Poincare sphere. 

(a) Points O, T and R describe the polarizations of o-, to- 
and ro-waves at the boundary of the twisted layer. The Mau¬ 
guin cone contains the point T on its axis, and the point 
O on its generator. Spherical triangle OTR has a right an¬ 
gle ZROT. The cathetus RO = H = 2^ lies on the “equator”. 
The cathetus TO = 0 = 2d is perpendicular to the “equator”. 
Acute angles are ZRTO = v and ZORT = po- The area is 
equal to spherical excess (OTR) = {Tr/2+v+po) — n = v — p. 

(b) Points E and Te describe polarizations of e- and fe-waves. 
The great circle GoG is perpendicular to the diameter TTe. 
The big circle TeGRT in the Gooch-Terry minimum coin¬ 
cides with the great circle TeGoOT. Distance to the Gooch- 
Terry minimum for the to-wave is determined by the phase 
u = ZGqG = D (TGoG), and for ro-waves it is determined by 
the phase p = D(GoGRO). 


Without changing the polarization the to-wave acquires 
a TT phase between the adjacent Gooch-Terry minima 
(see the section G 1 G 3 / 2 G 2 of the dispersion curve at 
Fig. [IKb)). The ro-wave gains the phase p (dispersion 
curve section G 1 RG 2 ) which is less than the phase v by 
the area of a spherical triangle O {OTR) = v — p. During 
the passage of the Gooch-Terry maximum this area has a 
critical value D {OTR) = 0 in accordance with Eq. (1241) . 

The cathetus RO=S = 2^ of the triangle OTR is given 
by Napier’s rules for right spherical triangles [s^ : 

tan S = — sin 0 tan t>. 

This equation determines the deflection angle ^ of the 
boundary ro-wave linear polarization from the LC direc¬ 
tor (Fig. [3]). FigureHKb) shows the arc TR rotating when 
V increases. The point R moves non-uniformly along the 
PS “equator”. The movement is faster near the Gooch- 
Terry maximum. In adiabatic approximation 0 <C 7 r /2 
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the point R jumps at the Gooch-Terry maximum so that 
^ Ri 0 for ro-wave and ^ r; 7 r/2 for re-wave. 

Exact correspondence between optical wave phases and 
characteristic space angles provides a visual support and 
a qualitative understanding of the phenomenon support¬ 
ing the validity of the obtained solution. 


F. Intermediate optical response presumption 

The motivation for this study was the debate about 
the direction of the transmission peak spectral shift in a 
TN layer. The intermediate optical response presumption 
was formulated for an anisotropic medium whose effective 
RI is between the ordinary and extraordinary RI, 

Uo < n < ne (31) 

Here are a few abstract examples supporting the inter¬ 
mediate optical response presumption. 

1. The average RI of nematic in the isotropic phase 

is [I!,[zi 

nio = (nl + 2n^) /3. 

2. In a homogeneous, uniaxially anisotropic medium, 
the extraordinary light wave propagating at an 
angle 9 to the optical axis has the following RI 
(Ref [3, Eq. (11.6)): 

rie (9) = noUe (0) / rig (0) cos^ 9 + sin^ 9. 

3. The thin sublayer of the TN layer has effective RI 
according to the normalization (I^Hl) . 

Non-additive response of a slab of sublayers leads to 
the geometric phase (1^ . As a result, the Mauguin for¬ 
mula (HI contradicts the intermediate optical response 
presumption 


nte,to =n±\j5n^ + Uto < Uo < Ue < Ute- 

The stated contradiction admits an experimental test. 
With an increase in the effective RI (i.e., increase in 
the optical length of the cavity), transmission peaks are 
shifted to the red. For the o-wave in untwisted struc¬ 
ture the effective RI is equal to the ordinary RI. Ac¬ 
cording to the intermediate optical response presumption 
(Eq. dsn) the twisted effective RI shifts towards the ex¬ 
traordinary RI; in other words, the effective RI increases. 
This predicts the redshift AAtss > 0 for spectral peaks. 
In contrast, the Mauguin formula predicts the blue shift 
AAtss < 0 (1271) . 


IV. EXPERIMENT 

A. LC orientation 

An experimental study of the shift of the o—polarized 
spectral transmission peaks to shorter wavelengths was 
carried out inside a LC-FPC. The Fabry-Perot cavity 
consisting of two dielectric mirrors (Fig. [6l bottom row) 
is treated as a LC cell. It was filled with the nematic LC 
4-methoxybenzylidene-4’-n-butylaniline (MBBA) doped 
with the cationic surfactant cethyltrimethyl ammonium 
bromide (CTAB) in the weight ratio 1: 0.003. The cavity 
gap is 7 pm. CTAB molecules within MBBA dissociate 
into Br-anions and surface active CTA-cations. The lat¬ 
ter are adsorbed at the surfaces of aligning layers and, 
at the sufficient concentration, can form a layer spec¬ 
ifying the homeotropic coupling condition for nematic 
molecules [s^.[^. 

A multilayer mirror coating consists of 6 layers of zir¬ 
conium dioxide (Zr02) with RI of 2.04 and a thickness 
of 55 nm and 5 layers of silicon dioxide (Si02) with RI 
of 1.45 and a thickness of 102 nm, alternately deposited 
on the surface of a quartz substrate. The alternating 
layers produce the reflection band in the range of 420- 
610 nm. Thin (~ 150 nm) ITO-electrodes were deposited 
on the npper layers of Zr 02 to apply the electric field nor¬ 
mally to the LC-FPC mirrors. Electrodes were covered 
with different polymeric alignment films by spin-coating 
to implement the initial homeoplanar director orienta¬ 
tion (Fig. HJa), bottom row). The top substrate was 
covered with a planar-orienting film of pure polyvinyl 
alcohol (PVA). The bottom substrate was covered with 
a PVA film doped with glycerin (Cl) compound in the 
weight ratio 1: 0.61. With the utilized CTAB concen¬ 
tration the surfactant ion molecules were adsorbed on 
the PVA-Gl film. The formed layer shields the planar- 
orienting effect of the polymer coating and provides the 
homeotropic anchoring for MBBA. Polymer films on both 
dielectric mirrors were rubbed unidirectionally to define 
the axis of easy orientation. The angle between the rub¬ 
bing directions at the top (RI) and bottom (R2) mirrors 
is 90°. 

Inside the LC-FPC cell at constant voltage value U 
= 4 V, the orientation transition from homeoplanar to 
twist director configuration is induced by ionic modifi¬ 
cation of surface adhesion [sl]. This transition is ac¬ 
companied by a modification in polarizing microscope 
optical texture of the LC-FPC cell under the parallel- 
and crossed-polarizer scheme (Fig.[BJ top row and middle 
row, respectively). For example, under crossed polariz¬ 
ers the optical texture of LC-FPC cells in the initial state 
is uniformly dark when the RI direction coincides with 
the transmission axis of the polarizer (Fig. UKa), middle 
row). The light transmission is increased when the volt¬ 
age is applied and the twist-configured FPC (TN-FPC) 
is formed (Fig. IHKb), middle row). However, under the 
parallel-polarizer scheme the TN-FPC transmittance is 
low (Fig. llKb), upper row) due to the rotation of the 
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(a) 









FIG. 6. Microphotographs of optical textures of LC-FPC under the parallel-polarizer scheme (top row) and crossed polarizers 
(middle row), and the respective diagrams of the LC director configurations (bottom row), (a) Homeoplanar orientation at U 
= 0 V; (b) twist orientation at U = 4 V. Polarizer (P) and analyzer (A) directions are represented by the double arrows. R1 
and R2 are rubbing directions of the top and bottom mirrors, respectively. 


polarization plane of linearly polarized light at an angle 
close to 90° after passing the TN layer. 

The experimental scheme excludes a signihcant im¬ 
pact of parasitic factors on the shift. The structure is 
twisted uniformly, because every sublayer has constant 
twist torque created by different rubbing orientations. 
Twist-structure is under the voltage of 4 V. However, 
this does not lead to any deformation of the structure. 
Firstly, in the volume of the cell the voltage is partially 
compensated by surface charges. Secondly, MBBA is ori¬ 
ented transversely to electric held. 


B. Cavity and its spectra 

Polarized transmission spectra of LC-FPC with two 
different conhgurations of MBBA director for normally 
incident light were measured using an Ocean Optics 
HR4000 spectrometer equipped with hber-optics (Fig.[71). 
The LC-FPC sample was placed inside the optical chan¬ 
nel with planar orientant (Ri || x) at the input mirror. 
A Clan prism (P) was used as a polarizing element with 
the linear polarization along the y-axis (i.e., orthogonal 
to the director on the input mirror). In this experimental 
scheme the o-polarized transmission spectrum was mea- 


























































13 



FIG. 7. Experimental setup for measurement of polarized 
transmittance spectra of LC-FPC. C - fiber collimator, GP - 
Gian prism, PD - photodetector. 

sured regardless of structural transformations in the ne¬ 
matic volume. Spectra were recorded at a fixed tempera¬ 
ture of 23.0° C. Thermal stabilization error was no more 
than ±0.2° C. The voltage was generated by a power 
supply (AKTAKOM ATN-1236). 


C. Comparison of experiment, simulation and 
analytic formulae 

Three measurements were averaged to calculate the 
experimental spectral shift (Fig. [8|). Experimental reso¬ 
lution of 0.25 nm was improved by fitting spectral peaks 
by Voigt contours. The confidence interval was calcu¬ 
lated as the corrected sample standard deviation mul¬ 
tiplied by the Student coefficient with 95% reliability: 
^0.05,2 = 4.30 2 7. 

The Berreman method [3^ [3^ was used for numerical 
simulations. The method evaluates polarization vectors 
and transfer matrices of dimension 4 to take into account 
the reflection between sublayers of the LC bulk. The TN 
layer was divided into 200 sublayers and the calculated 
spectral resolution was 0.01 nm. Some parameters were 
considerably tuned to match the experimental spectra. 
The thicknesses and RIs of amorphous layers constituting 
the dielectric mirror were taken for Si02: 83 nm and 
1.45; for Zr02: 66 nm and 2.02; for each ITO layer: 
117 nm and 1.88858 ± 0.022i taking the absorption into 
account; for the substrate RI: 1.45 and for the PVA RI: 
1.515; two thicknesses of PVA layers: 300 and 600 nm; for 
MBBA extraordinary RI: 1.737 and ordinary RI: 1.549, 
both with RI imaginary part: 0.00078i. The thickness of 
MBBA layer was 7980 nm; the twist angle was 80°. 

The material dispersion gave some minor changes in 
the spectra and the shift. The most notable changes were 
the absorption dispersion and the change in MBBA layer 
thickness by 40 nm. Therefore, the spectra illustrated 
in Figs. 8 and 9 were calculated assuming no material 
dispersion. 

Experimental transmission peaks are broadened no 



Wavelength, nm 



Wavelength, nm 

FIG. 8. Ordinary-polarized LG-FPC spectrum (a) and the 
spectral shift of transmission peaks (b). (•) - experimen¬ 

tal shift values, (o) - calculated shift values obtained by di¬ 
rect numerical simulation using the Berreman method, solid 
curve - analytical shift calculated from Eq. (1251) obtained by 
the Jones method, dashed curve - the shift AAxss calculated 
without mirror effect from Eq. (1271) . The parameters of the 
rightmost peak with A = 579.1 nm were taken to calculate 

Fig. SI 


more than 2 times larger than the calculated peaks 
(Fig. in]). This was achieved by additional flatness tun¬ 
ing of the cavity. Comparison of experimental spectra 
with simulations shows the transmission peaks shift to 
shorter wavelengths for the twisted configuration. The 
shift is less than the half-width at half maximum of a 
peak. Three peaks in the range of 480-500 nm have min¬ 
imal intensity and maximal shift dispersion. However, all 
14 points ht satisfactorily into the simulated dependence 
curve. 

Figure 8(b) shows by dots and circles the experimental 
and calculated values of the spectral shift, respectively. 
Analytical Eqs. (1^51) and (l?7l) are shown by the solid 
curve and the dotted curve, respectively. The reflection 
phase shift leads to simultaneous displacement of all the 
points along the solid curve. The latter analytic curve is 
affected by the reflection phase shift only near Gooch- 
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FIG. 9. Zoomed experimental (a) and calculated (b) trans¬ 
mittance spectral peaks for homeoplanar (blue) and twisted 
(magenta) configurations in LC-FPC. The twist leads to 
shorter wavelength shifts. 

Terry maxima. The experimental spectrum contains no 
Gooch-Terry maxima. So in analytical equations the re¬ 
flection anisotropy was ignored S^j, = 0. 

The analytic formula (051) slightly overestimates the 
shift. There are several possible reasons for this distinc¬ 
tion. Firstly, the difference of PVA layer thicknesses on 
mirrors leads to some error. Secondly, reflection between 
sublayers of LC bulk produces the difference between the 
Berreman and the Jones spectral shifts. The effective 
Jones RI (ITT)) is commonly used for the TN LC where 
the helical pitch is much larger than the wavel ength . 
The formula (ITTl) has been generalized [H, |H, [13 Issi) for 
cholesteric LC where the helical pitch is of the same order 
with the wavelength. The wave vector for the Berreman 
method: 



It gives the effective Berreman RI: 


rig = ±Jn^ + rfl ± \J+ {2n6n + 5rfi)^. (32) 


The Berreman RI ((32l) and the Jones RI (1141) contrast is 
evident in the graph scale despite the reasonable approx¬ 
imation <C J <C cr. 


V. CONCLUDING REMARKS 

To the best of our knowledge, the Mauguin phase shift 
formula m has had only implicit and indirect confirma¬ 
tions in polarization measurements. Such measurements 
are, for example, the measurement of spectral positions 
of Gooch-Terry minima and the optical activity mea¬ 
surement of cholesteric LC [t^. As usual in commer¬ 
cial TN cells TSS does not exceed 1 nm in accordance 
with Eq. (HZD. TSS is really hard to ascertain exper¬ 
imentally. For example, in the experimental spectra of 
TN-FPC with twist-homeotropic electric switching, the 
o-polarized spectral shift, as reported in [^ . was much 
higher than TSS Eq. (I?71) . The crucial factor was re¬ 
orientation of the parietal LC sublayers. They were not 
reoriented homeotropically up to the cell breakdown volt¬ 
age. 

The reported shift in TN-FPC can be observed directly 
without any polarizers. And the required measurement 
accuracy is achieved due to multiple interference in the 
cavity. The main drawback of the presented scheme is 
that untwisted structure maintains a constant RI for the 
o-polarized light only. The e-polarized spectral shift may 
be measured in the experimental scheme with the twist- 
planar to homogeneous-planar transition, which could be 
achieved using photoalignment material with reversible 
intermolecular bonds [H, [^ . 

The avoided crossing spectral shift caused by the mir¬ 
ror reflection mode coupling, the above-mentioned RSS, 
was described analytically and experimentally in the uni¬ 
formly twisted structure. The analytical expression (1261) 
for RSS near the Gooch-Terry maximum looks much 
easier than the one published before [s^. A general¬ 
ized Mauguin-Poincare rolling cone method allowed us 
to solve the problem geometrically, independently of the 
Jones and Berreman matrix formalisms. 

The spectral shift of transmission peaks should not be 
confused with a frequency shift or a frequency conver¬ 
sion. In the stationary linear problem the light frequency 
is not converted. Also we claim that the twist spectral 
shift characterizes not entirely the cavity but namely the 
twisted layer itself. The cavity just facilitates the mea¬ 
surement in that the twisted layer does not generate any 
transmission peaks. So what is actually shifted while 
twisting the layer? Obviously, shifted is the eigenwave 
phase when going out of the twisted layer. Therefore, 
the effective refractive index varies. This optical response 
can be measured without a cavity. For example, a polar¬ 
ization grating has considerable sensitivity to a minute 
change in refractive index, which permits the experimen¬ 
tal confirmation of the described phenomenon. 

The intermediate optical response presumption was 
formulated to determine the extreme values of the op- 
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tical response of the complex medium. Despite the ap¬ 
parent evidence this presumption may be violated. A 
well-known example is a composite of a few optical me¬ 
dia with the inhomogeneity scale much smaller than the 
wavelength. Its resonant optical response may exceed 
maximal values for every component. This intermedi¬ 
ate optical response violation for composites is explained 
by the Clausius-Mossotti contribution of spatial bound¬ 
aries between the components [U . Another interme¬ 
diate optical response violation for a twisted anisotropic 
medium is shown to be due to the GP contribution. This 
tiny GP contribution should not be confused with the 
GP of the zeroth-order adiabatic approximation. The 
last one is responsible for 7r-phase polarization conflict in 
the TT-twisted LG cell m- A visual GP representation 
is speculated as the area covered by the polarization tra¬ 
jectory on PS. Contributions of TSS and RSS are related 
to areas of spherical rectangles and triangles. 

The revealed tiny spectral shift in TN-FPC exists in an 


arbitrary anisotropic chiral medium, not only TN layer. 
No doubt, it must be accounted for a vast class of twist- 
polarization devices in high-precision engineering appli¬ 
cations, such as multiplexers, 3D displays, optical tweez¬ 
ers, holographic data storage, and diffractive optics. 
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